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ABSTRACT 

We present a comprehensive analysis of the spatial, kinematic, and chemical proper- 
ties of stars and globular clusters (GCs) in the 'ordinary' elliptical galaxy NGC 4494 
using data from the Keck and Subaru telescopes. We derive galaxy surface bright- 
ness and colour profiles out to large galactocentric radii. We compare the latter to 
metallicities derived using the near-infrared Calcium Triplet. We obtain stellar kine- 
matics out to ^ 3.5 effective radii. The latter appear flattened or elongated beyond 
^1.8 effective radii in contrast to the relatively round photometric isophotes. In fact, 
NGC 4494 may be a flattened galaxy, possibly even an SO, seen at an inclination of 
~ 45 degrees. We publish a catalogue of 431 GC candidates brighter than iq — 24 
based on the photometry, of which 109 are confirmed spectroscopically and 54 have 
measured spectroscopic metallicities. We also report the discovery of 3 spectroscopi- 
cally confirmed ultra-compact dwarfs around NGC 4494 with measured metallicities 
of —0.4 < [Fe/H] < —0.3. Based on their properties, we conclude that they are simply 
bright GCs. The metal-poor globular clusters are found to be rotating with similar 
amplitude as the galaxy stars, while the metal-rich globular clusters show marginal 
rotation. We supplement our analysis with available literature data and results. Using 
model predictions of galaxy formation, and a suite of merger simulations, we find that 
many of the observational properties of NGC 4494 may be explained by formation 
in a relatively recent gas-rich major merger. Complete studies of individual galaxies 
incorporating a range of observational avenues and methods such as the one presented 
here will be an invaluable tool for constraining the fine details of galaxy formation 
models, especially at large galactocentric radii. 

Key words: galaxies: haloes - galaxies: abundances, kinematics - galaxies: individual; 
(NGC 4494) 



1 INTRODUCTION 

Despite sustained efforts both on the observational and the- 
oretical fronts, an accepted self-consistent picture of galaxy 
formation and evolution has not yet emerged. Key processes 
involved during galaxy formation and evolution include (but 
are not limited to) galaxy mergers, the accretion and/or dis- 
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sipation of gas (whether early, via monolithic collapse, or 
merger induced), feedback processes such as stellar winds, 
supernova (SN) feedback or active galactic nuclei (AGN), 
reionisation, etc. 

Observational clues to the importance of these processes 
are crucial to our understanding of galaxy formation. For ex- 
ample, the relative importance of gas dissipation and energy 
feedback from stellar winds, SN and AGN as star formation 
quenching mechanisms needs to be constrained observation- 
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ally. In principle, this can be studied throu gh looking at ra- 



dial abundance gradients in a galaxy (e.R., Bekki fc Shioval 



1 19991 : iKobavashi fc Arimotdll^ : iHopkins et al.ll2009h . 

Also, the properties of progenitor galaxies involved in a 
merger can be probed by looking for distinct kinematic and 
morphological signatures in the merger rem nant galaxy (e.g., 
iLopez-Saniuan et al.l[201ol : lHopkins et al.i r2009. 20l d). espe- 
cially at large galactocentric radii. Indeed, Hoffman et al.l 
predict that the less relaxed intermediate and outer 
parts of a gas-rich (or wet) disk-disk merger remnant may 
retain the dynamical signature of the cold disk stars and 
original halo stars, respectively. In this model, a transition 
in the kinematic properties of the remnant is expected for 
1 < Te < 3 where is the effective radius. Moreover, as a 
complement to looking for these kinematic signatures, dry 
merging also leaves an imprint on the stellar populations 
of the remnant such t hat it washe s -out or weakens metal- 
licity gradients (e .g., IWhitd Il980l : iDi Matteo et al.l l2009l : 
iPipino et aklboiol ). 

Thus, in order to obtain a complete picture and under- 
standing of galaxy formation it is necessary to spectroscop- 
ically probe both the global kinematics and stellar popu- 
lations present in galaxies out to large galactocentric radii 
(i.e., > Ire). Unfortunately, the low surface brightness of 
galaxies in the outskirts hinders spectroscopic studies at 
large galactocentric radii. Spatially resolved kinematic and 
stellar population studies of the stellar light of large early- 
type galaxie s (ETGs) are typically confined to the inner 
< Ire (e.g 



> (-bl LrS) are typically conhned to tne inner 
Reda et al.1 l2007l : lEmsellem erahl |2007| . I2OIII : 



iKuntschner et al.ll201ol v! thereby probing less than half the 
stellar light and potentially missing the impor tant transition 
region previously discussed. To remedy this, IProctor et al.l 
l|2009l . hereafter P09) have developed a technique that takes 
advantage of the large field-of-view of the DEIMOS spectro- 
graph on Keck to extract galaxy light spectra from back- 
ground spectra (sky -I- galaxy light) of GC multi-slit obser- 
vation and derive spatially resolved ki nematics out to 3re 
Hoi) and metaUicities out to ~ 1.5re [Foster et al] l|2009l . 
hereafter F09). 

A good probe of early star formation are globular clus- 
ters (GCs). ETGs typically present a larger specific fre- 
quency (i.e., number per unit mass) of globular clusters 
(GCs) than late-type galaxies (e.g., Ivan den Berghl 1 19821 : 
I Ashman fc Zep3ll992l ). GCs are generally measured to be 
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Proctor et al.ll2008l) 



and thus probe the very earliest star formation episodes 
of their host galaxy. This known old age significantly di- 
minishes th e difficulties as sociated with the age-metallicity 
degeneracy llWorthevlliggil ) wherein old metal-poor popula- 
tions share certain photometric and spectroscopic properties 
with those of young metal-rich populations. Moreover, GCs 
are well approximated as single-stellar populations, which 
greatly simplifies their stellar population analysis and inter- 
pretation. 

Understanding GC formation is crucial to under- 
standing galaxy formation, especially the early epochs 
l|West et al.1 l2004l l. For example, the nearly ubiquitous 
colour bimodality of th e GC systems of most galaxies (see 
iBrodie fc Straderll2006l . for a review) is usually interpreted 
as a metallicity bimodality with the blue and red GCs be- 



ing metal-poor and -rich, respectively. This is a stringent 
constraint for galaxy formation scenarios that must pro- 
vide at least two formation episod es or me chanisms to ex- 
plain the colour bimoda lity (e.g., iLee et al. 2010; but see 
iMuratov fc Gnedinllioidh . 

A recent and extensive GC study is that of the kine- 
matics and stellar populations of a statistically sizeable frac- 
tio n of the GC sys t em of the nearby giant ETC NGC 5128 
by IWoodlev etal] l|2010l . a,b). They find that the major- 
ity of the GCs in both subpopulations are old with a sig- 
nificant population of young metal-rich GCs forming later. 
They conclude that these young metal-rich GCs may have 
been formed in a more recent merging event than the bulk 
of the GCs. Moreover, both metal-rich and -poor GC sub- 
populations are found to be pressure supported with only 
mild rotation for the metal-rich GCs. From this, they are 
able to infer that the GC system of NGC 5128 is consistent 
with a hiera rchical format i on in a sc enario similar to tha t 
proposed bv lBeaslev et al.1 l|2002l ) and lStrader et al.1 (|2005h . 
This study demonstrates the power of combined spectro- 
scopic and photometric studies of large samples of GCs to 
constrain the formation and assembly history of individual 
galaxies. 

While the most massive nearby elliptical galaxies are 
in principle easier to study, their formation may have been 
atypical due to their often special location at the centre 
of large potential wells such as galaxy groups or clusters. 
Therefore, in order to obtain a complete view of galaxy for- 
mation, it is necessary to avoid "special" environments and 
aim for more typical galaxy masses (i.e., ^ M*). To this end, 
we focus on the galaxy NGC 4494, which i s often described 
as an "ordinary elliptical" ga laxy (e.g., ICapaccioli et al.) 
1 19921 : iLackner fc Ostrikeij [2OI0I ) mainly based on its typi- 
cal light profile. Another aspect that makes it fairly ordi- 
nary, or average, include its intermediate density environ- 
ment. In deed, it has been described in the literature as either 
isolated (iLackner fc Ostrikedl201ol) or loose group member 
l|Forbes et al.|[l996l : Larsen et al.ll200ll ) as it is located at the 
edge of the Coma I cloud. It is neither a large nor a small 
galaxy with a stellar mass of ~ 10" M0 (see Table [T|. Its 
elliptical morphology is also typical as early-type galaxies 
(ETGs) may contain o ver 50 per cent o f the total stellar mass 
in the local Universe (iBeU et al.ll2003l). It c ontains an inner 
dust ring (r < 4 arcsec, Lauer et al.l 12005^. is quite round 
with an axis ratio of g = 0.87 (see Table[T| and shows a very 
smooth luminosity profile with a central cusp (|Lauer et al.l 
I2OO7I). Peculia rities include a kinematicall y decoupled core 
jBender et al.l J 994) or double maxima (jKrainovic et al.l 
I2OIII ) in the inner ~ 19 arcsec beyond which sustained ro- 
tation (Kot ~ 60 km s~^) is observed out to ~ 3re with a 
possible "pinchi ng" o r "fiattening" of the kinematics start- 
ing at ~ 1.5re l|P09r ). Other notable properties are a two 
order of magnitud es deficiency in X-ray lumin osity for its 
optical luminosity (|0'Sullivan fc PonmanI 2004) and a pos- 
sible deficiency in dar k matter ( Romanowskv et "all l2003l : 
iNapolitano et al.ll2009l . hereafter N09). 

Our approach is to study NGC 4494 in great detail to 
constrain the formation of "ordinary" elliptical galaxies. To 
this end we study its structure, kinematics and stellar popu- 
lations, as well as its GCs, using imaging and spectroscopy. 
We also include literature data. Our study is unique and one 
of the most complete studies of an individual "ordinary" el- 
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Figure 1. DSS image of NGC 4494 sliowing the positions of 
our science spectra. Circles and squares represent kinematically 
confirmed GCs and UCDs, respectively, while crosses show the 
position of galaxy light spectra. Green and black small symbols 
are used when the spectrum returned both CaT and kinematic 
information, while red symbols are used if only kinematic infor- 
mation could be extracted. Large black ellipses represent 3, 6 and 
9 re with the galaxy's global photometric position angle and axis 
ratio. The HST field-of-view is shown in white. A colour version 
of this figure is available in the online version of this article. 



liptical galaxy and its GC system to date. We present ar- 
guably the first large catalogue of GC recession velocities in 
an ordinary ETG. 

The paper is divided as follows: Section [2] presents the 
photometric and spectroscopic data. The surface bright- 
ness/density profile, kinematics, colours and metallicity dis- 
tribution of the stars and GCs are found in Section O These 
results and their implication for the key science questions 
described above are discussed in Section 2) Finally, we give 
a brief summary and our conclusions in Section [S] 



2 DATA 

2.1 Imaging acquisition and reduction 

Subaru Suprime-Cam (jMivazaki et al.l |2002| ) imaging of 
NGC 4494 is analysed to understand the GC system and 
light profiles of NGC 4494. We obtained a g-band observa- 
tion of NGC 4494 during a Gemini time exchange program 
(GN-2008A-C-12) on the night of 2008 April 2. Two addi- 
tional bands (r and i-bands) were acquired on a later date, 
2010 April 4, through Keck time exchange. The total ex- 
posure times are 805, 365, and 540 seconds in the g, r and 
i-bands, respectively. The seeing conditions in the respec- 
tive bands are: 0.63, 0.56, and 0.58 arcsecs. Suprime-Cam 
data are prepared for analysis using standard imaging reduc- 
tion tech niques and a modified version of th e sdfred data 
pipeline (|Yagi et al.ll2002l : lOuchi et al.ll2004h . The Suprime- 



Cam field of view covers a ~ 36 x 29 arcmin region centered 
on NGC 4494. 

Suprime-Cam photo metry is boots trapped to the Sloan 
Digital Sky Survey DR7 (|Abazaiian et al. 2009) photomet- 
ric system using point sources with 19 < i < 21.5 mag- 
nitudes. The estimated g-, r- and i-band systematic uncer- 
tainty due to this calibration is 0.004, 0.005 and 0.004 mags, 
respectively. All photometry is Galactic extinction-corrected 
according to ISchlegel et al.l l|1998il . 

The images are prepared for GC analysis by first mod- 
elling with IRAF/Ellipse and subtracting the galaxy light 
profiles. Ellipse is set to allow the center, position angle 
and ellipticity to vary. A bright blue star is ~ 6 arcmin away 
from NGC 4494 to the NNE. The scattered hght from this 
star extends to a radius of ~ 3.5 arcmin, thus IRAF/Ellipse 
is also used to model and subtract its light in the g, r and 
i-bands. Because the star and galaxy profiles overlap signif- 
icantly. Ellipse is performed and subtracted iteratively on 
the light profiles 3 times. The final image products show a 
constant background value across the field. 

A catalogue of GC candidates is constructed from the 3 
Suprime-Cam mosaics. At the distance of NGC 4494, GCs 
are unresolved and appear as point-sources on the images. 
IRAF/DAOPHDT/FIND is used to locate objects on the field 
deviating by 3.8a from a global background level. The dif- 
ference between two different aperture magnitudes is used 
as a way to identify a nd remove extended sources from the 
GC catalogue (see e.g. ISpitler et al.lboOSi '). This is done for 
detected objects in each image. Aperture corrections are ap- 
plied to the point sources and the photometric zeropoints 
derived from the above procedure are used. The three sepa- 
rate photometry catalogues are combined using a matching 
threshold of 1 arcsec. 

To supplement the Suprime-Cam catalogue, an existing 
Ifubble Space Telescope (IfST) Wide-field Planetary Cam- 
era 2 (WFPC2) GC catalogue is incorporated into the GC 
analysis. The WFP C2 catalogue is pro vided by S. Larsen 
and is described in iLarsen et al.l (|200ll '). See Section [S^T] 
for details on our selection of GC candidates. 

2.2 Spectroscopy acquisition and reduction 

Spectra were obtained as part of the SMEAGOL survej0 
on two separate dark nights: 2008 April 8 (hereafter Night 
1) and 2009 March 23 (hereafter Night 2) using the DEep 
Imaging Multi-object Spectrograph (DEIMOS) on the Keck 
telescope. Three slit-masks were observed during Night 1 
under good seeing conditions (FWHM ~ 0.7 arcsec) and 2 
slit-masks were observed during Night 2, this time under 
variable seeing conditions (0.8 < FWHM < 1.3 arcsec). The 
gala xy light data only for Night 1 were pres ented in both 
|P09l and lF09l . In addition to the data used in |P09l and lFOOl . 
the final dataset used here includes roughly twice as many 
galaxy light spectra and the whole set of GC spectra. The 
1200 1 mm~ grating centered on 7800 A was used together 
with 1 arcsec slit width. This setup yields a resolution of 
A A ~ 1.5 A and allows for the coverage of the Calcium 
II Triplet (CaT) spectral region (~8400-8900 A). A total 
of 3 and 4 half-hour exposures were taken yielding a total 

^ http:/ /sages. ucolick.org/surveys. html 
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M49 


E2 


163 


0.81 


15.1 ± 0.7 


104 


-21.5 ±0.16 


-25.5 ± 0.1 


3.1 ± 0.3 


997 ±7 


294 ± 3 


M60 


E2 


108 


0.81 


16 ± 1 


69 


-21.2 ± 0.2 


-25.2 ± 0.2 


2.4 ± 0.4 


1117± 6 


335 ±4 


M87 


cD 


152 


0.86 


15 ± 1 


95 


-21.3 ± 0.1 


-25.1 ± 0.1 


2.2 ± 0.2 


1307 ± 7 


335 ± 5 


NGC 1399 


El 


150 


1.00 


19 ± 1 


80 


-20.8 ± 0.3 


-25.0 ± 0.2 


2.0 ± 0.4 


1425 ± 4 


342 ±6 


NGC 1407 


EO 


60 


0.95 


27 ± 3 


70 


-21.4 ± 0.3 


-25.4 ± 0.3 


3.0 ± 0.9 


1779 ± 9 


272 ±6 


NGC 3379 


El 


68 


0.85 


9.8 ±0.5 


35 


-19.7± 0.1 


-23.7 ± 0.1 


0.59 ± 0.06 


911 ± 2 


209 ± 2 


NGC 4494 


El-2 


173 


0.87 


15.8 ±0.9 


49 


-20.4 ± 0.2 


-24.2 ± 0.2 


1.0 ± 0.2 


1344 ± 11 


150 ±4 


NGC 4636 


EO-1 


143 


0.84 


13.6 ±0.9 


89 


-20.2 ± 0.2 


-24.2 ± 0.1 


0.95 ± 0.08 


938 ±4 


203 ±4 


NGC 5128 


SO 


43 


0.89 


3.9 ±0.3 


305 


-20.2 ± 0.2 


-24.0 ± 0.1 


0.78 ± 0.07 


547 ±5 


120 ± 7 



Table 1. Properties of nine massive ETGs with significant GC kinematic samples from the literature. Hubble types (column 2) are as per 
NASA/IPAC Extragalactic Database (N ED). Position ang les and axis ratios (columns 3, 4) are from 2M ASS. Distances (colu mn 5) are 
based on surface brightness fluctuations llTonrv et al.ll2001 ) and include the distance moduli correction oflje nscn ct al.l ||2003| '). Effective 
radii (column 6) are taken from the Third R e ferenc e Catalogue of Bright Galaxies (RC3, Ide Vaucouleurs c t al. 199l]) for all galaxies 
except NGC 5128, which is from lDufour et al] l|l979h . B- and ii'-band absolute magnitudes (columns 7, 8) are calculated from RC3 and 
2MASS apparent magnitudes, respectively, and using the distances quoted in c olumn 5. Stellar masses (column 9) are calculated from 
the K-hand magnitude of column 8 assuming a M / Lk ratio corresponding to the lVOsI SS P of age 10 Gyrs and solar metallicity. Systemic 
velocities (column 10) are from NED. Central velocity dispersions (column 11) are as per lPaturel et al. I 1I2OO3I '). 
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Figure 2. Example GC spectra (black) and fitted pPXF tem- 
plates (red) for a bright (top panel), typical (middle panel) and 
faint (lower panel) GC in our sample. Shaded wavelength bands 
are regions affected by significant skyline residuals. GC apparent 
io magnitudes are given. 
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Figure 3. Example galaxy light spectra (black) and fitted pPXF 
templates (red) for a high (top panel), typical (middle panel) and 
low (lower panel) signal-to-noise spectrum. Shaded wavelength 
bands are regions affected by significant skyline residuals. Galaxy 
iQ surface brightness is given. 



exposure time of 1.5 and 2 hours per mask for Night 1 and 
2, respectively. Fig. [T] shows the positions of the slits that 
returned useful science spectra. 

The DEIMOS data are reduced using the idl spec2d 
data reduction pipeline provided online. Flat-fielding using 



internal flats, wavelength calibration using ArKrNeXe arc 
lamps, as well as the local sky subtraction are performed 
within the pipeline. The pipeline outputs the GC spectra 
with their corresponding fully propagated variance arrays 
as well as the subtracted background or 'sky' spectra for 
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each slit. Fig. [2] shows example GC spectra for a range of 
signal-to-noise ratios. 

We use the Stellar Kinematics with Multiple Slits 
(SKiMS) technique described in £09 and FOft to extract the 
galaxy light spectra from these background spectra. Indeed, 
a background spectrum is essentially the sum of the pure sky 
and the galaxy light spectra. Therefore, if one has a good 
estimate of the sky's contribution to the background spec- 
trum, it is possible to extract the galaxy light spectrum. For 
each mask, we obtain a high signal-to-noise estimate of the 
'pure' sky spectrum by summing and normalising several 
carefully selected background spectra at large radii (typi- 
cally -^6-7 re) where residual galaxy light is insignificant. 
This normalised sky spectrum is then scaled for each back- 
ground spectrum using the sky scaling index defined in lFOif I 
as the excess flux in the central passband with respect to the 
continuum level. The central passband of the sky scaling in- 
dex (8605.0-8695.5 A) is measured between the blue (8478.0- 
8489.0 A) and red (8813.0-8822.0 A) continuum passbands. 
This scaled sky spectrum is subtracted from the background 
sp ectru m to obtain a galaxy light spectrum. As mentioned 
in IFOQI . this method yields a flnal continuum level accurate 
to 0.7 per cent of the noise in the skyline residuals in the sky 
index definition region. The majority (~ 90 per cent) of the 
amplitude of these skyline residuals are likely due to the in- 
trinsic difficulties associated with non-local sky subtraction 
techniques such as variations of the sky spectrum across the 
DEIMOS field-of-view and with time. The additional ~ 10 
per cent of the amplitude of the skyline residuals in our sci- 
ence galaxy spectra are caused by small variations in the 
wavelength solution across the mask. Fig. [3] shows example 
galaxy light spectra for a range of signal-to-noise ratios. 



3 ANALYSIS AND RESULTS 

In this section we describe how we extract the photometric, 
kinematic and stellar population information for both the 
galaxy light and the GC system. We also give an overview 
of the method used to fit the kinematics of NGC 4494. These 
results are interpreted and discussed in Section U) 

3.1 Galaxy light 

3.1.1 Stellar light profile 

We first quantify the stellar light distribution of NGC 4494. 
Figure [3] shows the surface brightness profile extracted as 
described in Section [2. II in the g- and i-bands together with 
literature values. We fit Sersic profiles l|Sersid[T963l ') to the 
surface brightness profile {fJ,{r)) wit h geo metric radius r > 5 
arcsec to avoid the inner disk as per lNOOl . In practice, we fit; 

wher e 6„ = 1.9992n - 0.3271 (i.e. eq. 6 of ICraham fc Driver! 
I2OO5I ) for the Sersic index (n), the effective radius (r^) and 

^ The sky scaling index was redefined in IF09| in order to be tter 
avoi d spectral features and differs sligiitly from that of |P09| (see 
iFOgl . for details). 
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Figure 4. Surface brightness profiles for NGC 4494 as a function 
of intermediate radius (i.e., \/ab where a and b are the semi-major 
and minor axes, respectively). Filled circles are a comp ilation of 
our Subaru profile s and those f rom [ Lauer et al.l l|2005l l. Offsets 
were applied to the lLauer et al. ] ll2005f r V and / profiles to match 
our g and i profiles, respectively. Sersic fits to the outer regions 
(r > 5 arcsec) are shown as red solid lines. 



the surface brightness at the effective radius (/ie). The re- 
sulting fits are shown in Fig. [4] and Table [2] We measure 
values between 48 < re < 55 arcsec depending on the pho- 
tometric filter used. Throughout this work we use the liter- 
ature value of re = 49 arcsec ~ 3.76 kpc (Table [ij as it lies 
within the range of our measured values and is thus a good 
compromise. 

3.1.2 Stellar kinematics 

The stellar kine matic s of NGC 4494 have been probed out 
to large radii in |P09l using the SKiMS method. Here our 
galaxy spectral sample is roughly twice that of |P09l . We 
measure the velocity moments (recession velocity, velocity 
dispersion and Gauss-Hermite coefficients and /14) for all 
th e galaxy light spectra u s ing th e pPXF routine described 
in lCappellari fc EmsellernI l|2004l '). The pPXF routine uses a 
set of 13 stellar templates to determine the best overall kine- 
matic parameters and weighted combination of templates 
that minimises the residuals between the spectrum and the 
final resulting fit (see lCa ppellari fc EmscUcm 2004, for more 
details). The templates were observed using DEIMOS on the 
nights of the 2007 November 12-14 with a comparable instru- 
mental setup. They cover a range of stellar sizes (11 giant 
and 2 dwarf stars) and span spectral types from F to early 
M, which dominate around this spectral region. The fitted 
spectral range is limited to 8450-8750 A for stability and 
regions heavily contaminated by skylines are not fitted (see 
Fig- El - Each fit is carefully inspected for quality control. Un- 
certainties on the measured velocity moments are estimated 
using Monte Carlo methods. For each spectrum, we ran- 
domly reshuffle the residuals between the best fit template 
and the original spectra in the wavelength region fitted by 
pPXF before re-fitting. This is repeated 100 times for each 
individual spectrum. We use the standard deviation on the 
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Sample 


Us 1 Ne 


■r, 1 Re 


n 






(mag) 1 (#arcmin~^) 


(arcsec) 




(#arcmin 


galaxy g 


22.81 ± 0.03 


54.6 ±0.8 


3.72 ±0.09 




galaxy r 


21.90 ±0.03 


53.0 ± 0.6 


3.52 ±0.09 




galaxy i 


21.51 ± 0.02 


48.2 ±0.5 


3.47 ± 0.06 




all GCs 


5± 1 


100 ± 10 


1.7± 0.5 


0.28 ±0.02 


blue GCs 


2.2 ± 0.8 


138 ± 24 


1.8 ±0.7 


0.22 ±0.03 


red GCs 


3.1 ± 0.4 


83 ±6 


0.8 ±0.2 


0.053 ± 0.005 



Table 2. Results of Sersic fits to the the galaxy surface brightness profile as a function of geometric radius (Eq.[T] fitted for /Xe, fe and 
n) for points with r > 5 arcsec. Sersic fits to the GC density profiles (Eq. [6] fitted for N^, Re, n and bg) are also listed. 
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Figure 5. Variation in CaTpog and velocity moments of the 
galaxy light as a function of PA for spectra with r < Ire (black), 
Ire < r < 2re (red), and r > 2re (gray). Kinemetry fits are shown 
with solid lines for the observed recession velocity (Vobs) and ve- 
locity dispersion in each radial bin, as well as for /13 for spectra 
within Ire. Dashed lines represent the CaTpog saturation level 
(top panel), systemic velocity (second panel) and /13 = (bottom 
panel). The photometric position angle is PAp^ot = 173 degrees. 
This figure is available in colour in the online version. 



velocity moments for the 100 Monte Carlo realisations as 
our estimate of the random uncertainty. 

Fig. [5] shows the velocity moments as a function of 
position angle (PA) for our sample of spectra. Individual 
values can be found in Table lAll The galaxy stars show 
clear major-axis rotation. The amplitude of the rotation is 
roughly constant all the way out to > 2re and shows a flat- 
tening (i.e., axis ratio gkin = 0.42 ±0.06 is low) at large radii 
(r > 2re). The velocity dispersion decreases with radius. 
For the inner measurements, we find a clear trend with 
position angle although such a trend is not clearly visible 
beyond Ive- The fourth velocity moment (/i4) has constant 
amplitude at all radii within the uncertainties. 

In order to better understand the kinematic structure of 
NGC 4494, we use kinemetry to fit the model of an isotropic 
rotating ellipsoid (or inclined disk). Kinemetry is an analog 
to photometry where instead of fitting a model of the stellar 
light distribution , we f it a model of the kinematics (see e.g. 
fKrainovic et al ]|200i;|P0i). In what follows , we generalise 
and improve on the technique presented in IPOa to obtain 
kine metry fits to discrete and semi-discrete 2-dimensional 
data. iPOg have shown that kinemetry using this model pro- 
yides a good fit to th e inner regions probe d by SAURON 
lEmsellem et al.ll2004l ). In contrast to |P09l where the data 
were binned in circular annuli, here the data are binned 
in overlapping 'radial' (semi-major axis) intervals that, to- 
gether with the position angle and axis ratio, define 'ellipti- 
cal annuli '. The j^^ annulus (or bin) contains a subset of Nj 
observed recession velocity data points {Vobs,i)- We initially 
assume that the kinematic position angle and axis ratio of 
the galaxy coincide with the photometric values. The po- 
sition angle and axis ratio of the j*'^ elliptical annulus are 
iteratively refined to match the kinematics as the algorithm 
fits for the kinematic position angle (P^kin.j), axis ratio 
(gkin.j) and the amplitude of the rotation (Kot,j) using 
minimisation. 

In practice, we first define an inner elliptical annulus 
{j — 1) whose short semi-major axis length (oj) is the dis- 
tance to the closest data point and whose long semi-major 
axis length is aj + Aaj. We perform a minimisation of 
the data points contained within this elliptical armulus to 
the model. For the j"^ annulus, the is computed using 
the following equation: 



(Vo, 



Vr 



mod, i,j ) 



(2) 



where 
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mod,i ,j 



Vrot ,j 



1 + 



(3) 



and the ambivalent sign is positive (negative) if {PAi — 
f'^kin.i) is in the first or fourth (second or third) quadrants. 
In Equations[2]and[3l PAi and AV^^^g j are the position angle 
and the uncertainty on the recession velocity measurement 
of the i"' data point, respectively. The systemic velocity of 
NGC 4494 {Vsys = 1338.5km s"^) is measured by fitting the 
inner, highest signal-to-noise, data only prior to fitting the 
full dataset. 

For each iteration within the j*'^ annulus, a new set of 



degrees) for points beyond ~ Ire (i.e., red and grey points), 
suggesting that the kinematic flattening (i.e., low gkin) is 
present in the velocity dispersion also. While this dip is sug- 
gestive, we are unable to reliably fit the axis ratio of the 
velocity dispersion explicitly with the current dataset and 
assume that of the recession velocity moment. 

In practice, for the j^^ bin we fit the velocity dispersion 
(uj ) using x^-minimization where 



2 



= E 



fobs 



Act' 



(4) 



Eq. 2] is different in form from Eq. [5] because the velocity 



parameters {Vrot,j, PA^n,, , <7kin., ) is found and the values of dispersion is an even moment. In Eq.|l aobs,> and Aa', 



PA]iin,j and (jkin,j are used to update the shape and orienta- 
tion of the elliptical annulus to be used in the next iteration 
{aj and Aaj remaining fixed). Therefore, by construction 
this changes the selected member data points in the j^^ bin 
slightly for each iteration. After 15 iterations, we move out 
to the next 'radial' bin (i.e. {j + 1)) by increasing the length 
of the short semi-major axis such that flj+i = flj -I- (Aa)/3. 
The factor of 1/3 used here is chosen to increase the num- 
ber of radial sub-samples or bins in order to smoothly define 
the kinemetry radial profile. The semi-major axis bin width 
(Aflj+i) incr eases as a function of a-,+i according to a de 
Vaucouleurs (|de Vaucouleurslll953l ) surface brightness pro- 
file (i.e. oc l/(aj+i)^/*) in order to compensate for the lower 
signal-to-noise ratio at large galactocentric distances. With 
this setup, there are 15 overlapping bins, each containing 
~ 25 data points. This process continues until the radial 
extent of the data is covered. 

In order to obtain a /dof ~ 1 at all radii and partic- 
ularly in the inner regions we add an uncertainty of 5 km 
s~^ in quadrature to the Monte Carlo uncertainties (i.e., 
AV;'b,,i = V(AT4bB,>)2 + (5 kms-i)2, see Table EU- While 
the relative amplitudes of the uncertainties are crucial in 
order to obtain stable fits, this monotonic increase of the 
estimated uncertainties has a negligible effect on our fitted 
values. Nevertheless, it is an indication that either (1) there 
are velocity substructures at small galactocentric radii such 
that the model of a simple rotating isotropic ellipsoid is not 
appropriate or (2) the uncertainties are slightly underesti- 
mated. Indeed, the latter is likely since our quoted random 
uncertainties are obtained using Monte Carlo methods on 
our spectra and not independent measurements. For exam- 
ple, due to the finite (non-zero) length of the slits on the 
mask, their position angle, the effects of seeing, etc, indepen- 
dent measurements of the kinematics for the same "position" 
around the galaxy using either another instrument, setup or 
mask could yield slightly different values. For this reason, 
there could reasonably be some unaccounted for systemat- 
ics of order ~ 5 km s~^. Finally, two data points deviating 
from the best fit at the > 3(t level were excluded. 

The kinemetry fits for the higher order velocity mo- 
ments (velocity dispersion and ho,) are done in parallel in 
each annulus and using the values of PAkin.j and gkin.j ob- 
tained from the recession velocity kinemetry fits. In principle 
the position angle and axis ratio of these higher order ve- 
locity moments need not be equal to those of the recession 
velocity, but our data do not constrain these two parame- 
ters properly. Indeed, the velocity dispersion data in Fig. [6] 
show a hint of a dip along the minor axis {PA ~ 83 and 263 



obSjZ 



> km s" 



)2 are the measured recession velocity 

I"" spectrum and its associated random uncertainty 
(Auobs.O with an additional uncertainty of 8 km s~^ added 
in quadrature in order to obtain a \^ /dof « 1 as for the fits 
to Vohs above, respectively. 

A similar approach is applied to fit the amplitude of 
the third (H-j,) moment. Measurement uncertainties on indi- 
vidual velocity moments are propagated for the kinemetry 
fits using Monte Carlo methods. For each individual data 
point we have a measured uncertainty. Assuming that these 
uncertainties are Gaussian, we have a known distribution of 
possible measurements for each data point. Thus, for each 
point we randomly select from that distribution. Each time 
this is done for the entire dataset and we get a new set of 
possible measurements. The kinemetry algorithm described 
above is then applied to the new dataset. This exercise is 
repeated 100 times. From the resulting range of fits we de- 
termine the 68 per cent confidence interval of the best fit 
model allowed by the data. 

We also use Monte Carlo methods to verify that kineme- 
try can be reliably applied to sparse data. We create a series 
of rotating ellipsoid models by varying the input rotational 
velocity Kot.in = 50, 80 and 100 km s~^, and kinematic axis 
ratio Qkin.in = 0.5, 0.8 and 1.0. For each model, we randomly 
sample 115 data points within ~ 3re and assign measure- 
ment errors consistent with our observed errors at that ra- 
dius. This is repeated 25 times for each of the (K-ot,in, (Jkin.in) 
combinations. We verify how well the input parameters are 
then recovered using kinemetry. In general, we find that the 
standard deviation of the output rotational velocity varies 
between crv^ot out — 5 and 7 km s~^, while cTq-^^^^^^^ ^ 0.1 
for all gkin.in- We find that the fits are less stable when both 
Kot.in and gkin.in are low. We do not find other "unexpected" 
systematics. We conclude that kinemetry can be applied to 
sparse data successfully provided that there is sufficient ro- 
tation. 

The results of the kinemetry fits to the data are shown 
in Fig. [B] In general, we find good agreement bet ween 
the long-slit and p lanetary nebulae results of IN09| and 
ICoccato et"ai] l|2009l ) and our data from the SKiMS method. 
We have also tried fitting the PNe data from INOgI using our 
method for discrete velocities (see Section IB. 2. 31) a nd find 
good general agreement with ICoccato et al. I l|2009l ) except 
around 2re r < 3re where the fits did not converge due 
to an appare ntly l ow rotation. The trends described above 
and found by[P09| are confirmed. We find sustained major- 
axis rotation at all radii consistent with both the results of 
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Figure 6. Results of our kinemetry fits (solid lines) as a function of the major axis equivalent radius (Ire = 49 arcsec ^ 3.76 kpc) with 
68% confidence intervals (hatched regions). Black, dark blue and red represent fits to the stellar light, the blue and red GCs, respectively. 
The top panel shows the value of the per d egr ee-of-freedom for the fits of the first velocity moment only (i.e. Vrot)- Light blue, green 
and orange filled symbols show results of iNOOl and lCoccato et ahl ||2009|) from long-slit spectra, planetary nebulae (PNe) and galaxy light 
photometry, respectively. This figure is available in colour in the online version. 



IN09| and ICoccato et al.l (|2009t ) . We do not find evidence for 
minor-axis rotation. There is a slight decrease in velocity 
dispersion that levels off with radius at a value of ~ 100km 
s~^ beyond ~ 2.3re slightly at odds with the velocity disper- 
sion estimated from planetary nebulae but consistent with 
that of the GCs (see Section I3.2.3|l . The amplitude of the 



third velocity moment levels off beyond ~ 1.8re at a value 
of i?3 0.03. 

lEmseUem et al.l (l2007t ) classified ETGs as slow- or 
fast-rotators according to a luminosity- and radius- 
weighted value of the angular momentum proxy Xr = 
Vrot{r) / \/Vrot{rY + o{rY, where slow- and fast-rotators 
correspond to Xr < 0.1 and 0.1, respectively. These clas- 
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sifications were made for radii < Ire, but |P09| demonstrated 
that the picture could change dramatically with more radi- 
ally extended data. In the case of NGC 4494, we find that 
a local Xh ~ 0.35 (i.e. fast rotation) at virtually all probed 
radii, modulo a sharp di p around O.lre due to the kinemat- 
ically decoupled core. In lEmsellem et al.1 (|201ll 'l. NGC 4494 
is also found to be a fast-rotator using a new classifica- 
tion method that takes into account the ellipticity of the 
isophotes. 

Although not as pronounced as originally found in |P09l . 
a slight kinematic twist is visible beyond ~ 2re such that 
the kinematic position angle is different from the photomet- 
ric position angle at large radii. We measure a radial change 
in the axis ratio such that the velocity distributio n bec omes 
fiattened (i.e., small gkin) as originally detected bv lP09l . This 
behaviour is at odds with the measured isophotes that re- 
main much rounder even at large radii with no sign of sig- 
nificant diskiness. We test the robustness of the fiattening of 
gkin at large radii using three methods. The first is the above 
described nominal method in which we allow gkinj to freely 
vary and iteratively define the bins accordingly. Secondly, 
similarly to the binning method used in |P09l . we let the fit- 
ted Qkin.i vary but fix the bin shapes to Qphot , i.e. we do not 
let the kinematic axis ratio define the axis ratio of the bins. 
While this yields an overall larger /dof value, lower values 
of gkin,j are still marginally preferred at large radii. Thirdly, 
we fix Qkin.j = (Jphot for both the binning and the fitting 
and for all bins (j). This yields substantially larger /dof 
values than in the (first) nominal case. We conclude that 
the fiattening of the kinematics is most robustly detected 
in the first method but is still marginally detected when 
the kinematic axis ratio is allowed to vary inside selection 
bins whose shape have fixed axis ratios. On the basis that 
the only self-consistent methods are the first and the third 
and given that the third method yields significantly higher 
X^/rfo/ values than the first, we conclude that the former is 
preferable. We also tried removing the apparent "outliers" 
located around PA ~ and 360 degrees with recession ve- 
locities below 1570 km s~^ to see if they were causing the 
fiattening and the results did not change significantly. We 
conclude that the kinematic fiattening is robust (i.e., gkin is 
indeed low). 

Fig. [7] shows the relationship between measured and 
{Voha — 'Kiys)/o", which is a measure of the fraction of rota- 
tional over pressure support. The anti-correlation between 
/i3 and {Vohs, — Vbyb) /ct is indicative of a disk kinematic struc- 
ture. This may be the source of the fiattened kinematics. 



3.1.3 Stellar colours and metallicities 

We investigate the radial colour gradient of the stars in 
NGC 4494. By doing this, we obtain a rough estimate of the 
radial metallicity gradient by assuming a fixed old GC-like 
age at all radii. In other words, we attribute any colour vari- 
ation to changes in the metallicity (i.e. ignoring age effects). 
We convert (g — i)o colours into approximate metallicity us- 
ing th e empirical linear relationship derived bv^ Sinnott et al.l 
||2010D for GCs in NGC 5128. This is shown in Fig. H The 
galaxy stellar light is generally red and shows an overall 
colour /metallicity gradient. Because all the visible du st is 
contained within the inner 4 arcsec iLauer et aklbOOST ). we 
will assume that dust is not affecting our inferences on the 
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Figure 7. Relationship between the higher order Gauss-Hermite 
moments hs and (Vobs — Vsys)/cr. Point sizes are proportional to 
the signal-to-noise ratio of the spectra. The anti-correlation be- 
tween hs and (Vobs — Kiys)/"" suggests a disk structure is present. 
Typical error bars shown in the upper left corner. 



metallicity at the radii probed here. Therefore, between 10 
and 70 arcsec (or ~ 0.2re < r < 1.4re), the colour pro- 
file suggests a moderate metallicity variation with radius of 
—0.17 ± 0.02 dex per dex. Assuming a younger age would 
increase the colour- inferred metallicity. Beyond 2re the pho- 
tometric uncertainties start to dominate. 

We use the CaT as a spectroscopic proxy of metallicity. 
The CaT has been employ ed as a metallicity indicator for re- 
solved red giant stars fe.g..lDiaz e t al.lll989l : [jorgensen et al.l 
1992; Koch ot al. 2006), i ntegrated light spectra of Galactic 
('e.g.. lBica fc Alloinlll987l:lArmandroff fc Zinnlll988l ') and ex- 
tragalactic ( Foster et al.l 20101. hereaft e r FIO ) GCs as well 
as galaxies (e.g., Cenarro et al.l [200I : |F09| ). We compute 
the CaT index valu e for each galaxy spectrum using the 
index definition from IFO^ (i.e. CaTpog). The CaTpog is cor- 
rected for velocity dispersion broadening and converted into 
\Fe/H] using the sing le stellar population (SSP ) mo dels of 
IVazdekis etahl (|2003l . hereafter V03) following iF09l . Once 
again, we emphasise that the inferred galaxy starlight metal- 
licities are uncertain due to the above caveats. Fig. [8] shows 
the CaT index as a function of the semi - majo r axis equiva- 
lent radius (as per equations 3 and 4 of|F09|). The CaTpog 
gradient is found to be essentially fiat at all probed radii in- 
dicating an u ndetectably small radial change in metallicity. 

The lVO j SSP models predict that the CaT features sat- 
urate at a metallicity of roughly [Fe/H] ~ —0.5 such that 
there is a maximum allowed CaTpod value of 6.18 43- This 
theoretical behaviour has not been confirmed observation- 
ally. Our galaxy CaTpog data are consistent with a satu- 
ration around 6.18 A. Because our data scatter about this 



This is different from the saturation value for CaTpio because 
CaTpo9 is an altogether different index. Therefore, the absolute 
values of the two CaT indices used in this work for GC and galaxy 
light spectra cannot be directly compared. 
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Figure 8. Panels (a) and (b) show the CaTpog index and (g — i)o colour gradients of the NGC 4494 stellar light, respectively. Thin 
dashed lines show the saturation limit p redicted by the | V 03lSSP models (i.e., CaTpog = 6.18A or [Fe/H] ~ —0.5 dex). Thick solid line 
in Panel (a) is a rolling average (see e.g., ISawilowskvll2007l'l using bins of 8 data points. Panel (c ) is a colour histogram of the photometric 
sample of GCs with io < 23.5 and {g — i) conversion to metallicity as per lSinnott et al] l|201Cll '). 



Passband 


CaTpiO 

(A) 


C aTpog 

(A) 


Gal 


8490.0-8506.0 


8483.0-8513.0 


Ga2 


8532.0-8552.0 


8527.0-8557.0 


Ga3 


8653.0-8671.0 


8647.0-8677.0 


Shared continuum 




8474.0-8483.0 
8514.0-8526.0 
8563.0-8577.0 



8619.0-8642.0 
8680.0-8705.0 



Table 3. GaT index definition for galaxy stellar light spectra 
(CaTpog) and GG spectra (CaTpio). 



limit we cannot confidently use thelvoal models to convert 
our CaTpog values into metallicity. It is possible however 
that, as the colours suggest, the metallicity variation in the 
radial range probed is small (only A[Fe/H] ~ 0.2 dex) and 
thus a variation in CaTpog may be hard to detect within 
the uncertainties. Indeed, between 20 and 80 arcsec, where 
the vast majority of our CaT data lies, both the colours and 
CaTpog values are consistent with the predicted saturation 
limit. 



3.2 Globular cluster system 

3.2.1 GC colours and candidate selection 

As described in Section [2TT] we have a combination of both 
HST/WFPC2 and Subaru/Suprime-Cam imaging available. 
Various selections are applied to identified point sources in 
order to avoid contamination by foreground stars and back- 
ground galaxies. We first apply an upper size cut in both the 



HST and the Subaru imagefl in order to remove clearly ex- 
tended objects such as background galaxies. We also apply 
a spatial cut removing all objects beyond 8 arcmin (~ lOr^) 
from the galactic centre to avoid further contamination that 
dominates beyond that radius (see l3.2.2ll . Our colour selec- 
tions are as follows: (1) whenever HST imaging is available, 
GC candidates are selected based on [V — I)o colours only. 
(2) If HST imaging is not available, the selection is based on 
a colour-colour cut in (r — i)o vs {g — r)o space from Subaru 
photometry (see Fig. |9]). The final master GC catalogue of 
Subaru and HST photometry includes 431 selected GC can- 
didates brighter than io — 24. This catalogue is available 
online (also see Table |4ll. 

Selecting bright objects in common between the HST 
and Subaru images beyond a galactocentric radius of 0.5 
arcmin to avoid Subaru image artefacts, we obtain the fol- 
lowing colour conversion: 

(5-i)o = 1.23 X (l/-7)o-0.25 (5) 

with an rms= 0.03 5 using a linear bisector fit 
(jPeigelson fc Babulll992l ). We use this conversion to convert 
HST {V — I)o colours into {g — i)Q to merge with the Subaru 
photometry. 

The final colour magnitude diagram (CMD) for all our 
GC candidates with io < 24 is shown in Fig. 1101 Three 
bright objects (io < 19.4 mag) have recession velocities con- 
sistent with that of NGC 4494 but their absolute luminos- 
ity (Afi < —10.9 mag) suggests they fit within the defini- 
tion of ultra-compact dwarfs (UCDs) , or equivalently, dwarf- 
globular transition objects (DGTOs, see Section [3.2.511 . All 
three UCDs are centrally located (within 2 arcmin or ~ 4re 
of NGC 4494's centre). Other bright objects (io < 19.4 mag) 

* The size cut based on Subaru is only applied to objects beyond 
galactocentric radii of 0.5 arcmin because of crowding issues. 
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Table 4. Catalog of photometrically selected GC candidates with io < 24. Columns 1 and 2 give the position in right ascencion and 
declination {J2000), respectively. Columns 3 to 7 are the Subaru/Suprime-Cam and HST/WFPC2 photometry. The full table is available 
in the online version. 



(hh:mm:ss) (hh:mm:ss) 
(1) (2) 



SO 
(mag) 

(3) 



ro 
(mag) 
(4) 



^0 
(mag) 

(5) 



Vb 
(mag) 

(6) 



(mag) 
(7) 



12:31:23.70 
12:31:07.44 
12:31:38.83 
12:31:08.79 
12:31:18.50 
12:31:36.43 



25:40:12.90 
25:40:15.96 
25:40:17.01 
25:40:27.99 
25:40:29.10 
25:40:34.57 



24.072±0.023 
24.362±0.026 
24.922±0.032 
23.805±0.022 
23.377±0.020 
22.326±0.019 



23.359±0.020 
23.672±0.022 
24.271±0.030 
23.237±0.018 
22.823±0.016 
21.864±0.015 



22.981±0.019 
23.319±0.021 
23.830±0.028 
22.991±0.020 
22.562±0.018 
21.662±0.016 
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Figure 9. Visual representation of our applied colour-colour se- 
lection for objects with galactocentric radii < 8 arcmin brighter 
than io = 23.5. Filled green star symbols are spectroscopically 
confirmed GCs while hollow red stars show our photometrically 
selected GC candidates. Solid line is the inferred GC sequence 
(i.e., (g — r)o = 1.450(r — i)o + 0.142). Spectroscopically con- 
firmed stars (orange asterix) and emission line galaxies (orange 
filled circles) observed either as contaminants or fillers are shown. 
This figure is available in colour in the online version. 



with colours consistent with the NGC 4494 GCs are uni- 
formly distributed across the Subaru/Suprime-Cam field-of- 
view (19x15 arcmin) suggesting that they are likely con- 
taminant foreground stars. For this reason we apply a bright 
magnitude cut at io = 19.4. 

We check for the usual colour bimodality within our 
GC candidates in the combined HST and S ubaru (q — i)o 
colou r distributions. We apply the KMM (Ashman et al. 
Il994h test for the heteroscedastic (unequal widths) case to 
our GC candidates with io < 23.5. The colour histogram 
and fit results are shown in Fig. IIUI The returned p-value 
is < 0.0001, suggesting that a bimodal colour distribution 
is strongly favoured over a unimodal one at the > 99.99 



per cent level. The peaks for the blue and red GCs are 
at (g - i)o = 0.84 and 1.07, with w idths of cr = 0.084 
and 0.051, respectively. For comparison. iLarsen et all (1200 ll ) 
find peaks at {V — I)q = 0.90 and 1.10, corresponding to 
(g — i)o = 0.86 and 1.10, for the blue and red subpopula- 
tions when using the homoscedastic test on the HST data 
only. These small discrepancies are likely due to the intrin- 
sic differences between the datasets (HST vs Subaru), pos- 
sible intrinsic radial colour gradients (e.g., Harris 2 0091) and 
metho ds (homoscedastic vs heteroscedastic). iLarsen et al.l 
(I 2OOII ) also infer a relatively equal number of GCs in each 
subpopulation. We choose to apply a nominal colour split at 
{g—i)o ~ 0.99, corresponding to (V^— 7)o ~ 1.01 based on the 
KMM results in order to sep arate the blue and red GCs. Us- 
ing the lSinnott et all (|201G| ) conversion from (g — i)o-colour 
to metallicity suggests that our colour split corresponds to 
a metallicity of [Fe/H] ~ —0.68 dex. 



3.2.2 GC spatial distribution 

We construct projected surface density profiles of the 
NGC 4494 GC system candidates. Independent photometric 
selection criteria (as described in Section r3.2.ip are applied 
to both the HST and Subaru datasets to generate GC cat- 
alogues for this analysis. Between 1 and 4 arcmin, there 
are reliable surface density data for both Subaru and HST 
GC catalogues. Within the uncertainties the Subaru profile 
agrees with the HST profile, as shown in Figure [TT] 

We fit the GC density profile (7V(r)) with a Sersic pro- 
file (fSersic 1963.) similar to that commonly done for galaxy 
surfac e density profile s and recently extended to GC sys- 
tems JPeng et al. 201ll). We fit t he following variation of eq. 
1 from lGraham fc Driver! l|2005l ): 



N{r) = Ne exp (-&„ [{r/Re}'^" - l)] + bg, (6) 

where 6„ — 1.9992n — 0.3271. Free parameters recovered are 
the Sersic index (n) of the GC system, which is a measure 
of the steepness of the profile, the effective radius of the 
GC system {Re), which gives us a measure of the extent of 
the GC system, the surface density at that radius {N^) and 
the background or contamination level (bg). The results of 
the fits are shown in Fig. [11] and Table (2] For comparison 
with previous studies, we also fit a Power-law profile (i.e., 
N(r) oc r''). We obtain slopes of a = -1.7 ± 0.2, -1.8 ± 0.2 
and — 2.2±0.3 for all, blue and red GCs, respectively. As with 
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Figure 10. CMDs (top panel) for our GC candidates beyond 0.5 
arcmin from the galactic centre selected based on Subaru (hol- 
low circles) and HST (filled circles) imaging. All spectroscopically 
confirmed GCs around NGC 4494 are shown as crosses. Typical 
Subaru error bars are shown on the right hand side of the top 
panel. Lower panel shows histograms for our candidates brighter 
than io = 23.5 (hollow) with KMM fits for each/the sum of the 
two subpopulations as solid/dashed grey line(s) and confirmed 
(filled histogram) GCs. The colour distribution of the confirmed 
GCs is representative of that of the candidates. Dotted line rep- 
resents our fiducial colour split. 
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Figure 11. Surface density profiles of the NGC 4494 GC sys- 
tem. All, red and blue subpopulations are shown in black, red 
and blue, respectively. Open and filled symbols show the HST 
and Subaru data, respectively. Sersic fits with and without back- 
ground/contamination constant (see text) are shown as dotted 
and solid lines, respectively. Also shown are scaled and offset 
Sersic profile fits to the galaxy surface brightness profiles in the 
g- and j-bands (^) as labelled (see Section 13.1.11 and Fig.|4]l. This 
figure is available in colour in the online version. 



20111: ISnvder et al.ll201ll ). Therefore, using equation 3 of 
Harris fc Harrisl ~ 201ll '). we estimate that the mass of the 
central black-hole in NGC 4494 is (1.6 ± 0.2) x 10* Mq. 



other gal axies (e.g.. iForbes et al.|[l99i : iDirsch et al.ll2003l : 
iBassino et al.. 2009 ). the red GCs are more centrally concen- 
trated than the blue GCs. This is also inferred from their re- 
spective effective radii (i?e,biuc = 138 arcsec and Reared ~ 83 
arcsec) . 

We use 500 Monte Carlo realisations to estimate the 
total number of GCs based on the best fit Sersic param- 
eters and their covariance matrix by separately integrat- 
ing the density profiles for all, the blue and the red GCs 
to infinity. We also explicitly account for completeness and 
contamination. This yields an estimated number of GCs of 
392 ± 49, 324 ± 74 and 125 ± 10 for all, red and blue GCs, 
respectively. Therefore, assuming my — 9.70 (RC3) as to- 
tal galaxy magnitude and using ot her variables as per Table 
T] th e specific GC frequenc y Sn ([Harris &: van den Berghl 
1983) is 1.2 ±0.3, while Tn (|Zepf fc Ashmanll 19931 ) is 4 ± 1. 
This is somewhat low for galaxies of s imilar stellar mass 
but still within t he observed range (see lSpitler et aLll2008l : 
iPeng et al.|[2008l ). 

It has recently been shown that the correla- 
tion between total number of GCs and the cen- 
tral black-hole mass (|Spitler fc Forbes! 120091 ') is tighter 
than other pr eviously observed correlat i ons with black- 
hole mass (see iBurkert fc Tremaind l2010l : [Harris fc HarrisI 



3.2.3 GC kinematics 

For the candidate GCs that have available spectra with suf- 
ficient signal-to-noise ratio, we first measure their recession 
velocity. We use the iraf task rv.fxcor to cross-correlate 
the spectra with the 13 stellar templates described in Sec- 
tion 13.1.21 Thus, for each GC we have 13 measured reces- 
sion velocities together with output uncertainties from the 
FXCOR routine. The recession velocities reported in Table 
IA2I are the average measured recession velocities using all 
13 template cross-correlation results. The quoted uncertain- 
ties for the GC recession velocities are the maximum be- 
tween 5 km s~^ or the average output uncertainties given 
by FXCOR, which correspond to the average width of the 
cross-correlation peaks, added in quadrature to the stan- 
dard deviation among the templates, which is an estimate 
of the systematics. Some of our spectra cover sufficiently 
blue wavelengths to use the Ha feature at 6563A to confirm 
the measured recession velocity. Some of our spectra were 
from Galactic stars and background emission line galaxies. 

^ We fix Pj4]jin = PApijot = 173 degrees to get a better handle 
on the red GC kinematics since the rotation is low (see text). 
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Table 5. Compilation of salient GC kinematic properties in the literature for selected large ETGs (column 1). Rotation amplitudes for 
all, blue and red GCs are given in columns 2, 5 and 8, respectively. The position angle of the rotation for all, blue and red GCs are 
given in columns 3, 6 and 9, respectively. Photometric position angles for each galaxy can be found in Table [l] The velocity dispersion 
for all, blue and red GCs are shown in columns 4, 7 and 10, respectively . The numb er of GCs in the kine matic sample for e ach st udy is 
shown in column 11. R eferenc es in column 12 correspori d to ICoie et al.l ||2003|. C03).lHwane et al.l ||200£ , H08 ) , ICotij et al.l ll200ll. COl), 
ISchuberth et al.l l|2010l . SIO), iRomanowskv et al.l l|2(309l . R09). iBergond et all l|2006l. B06). iLee et al.l j20ld. LIO) and IWoodlev et"al] 
||2010| . WlOa). 



Table 6. List of spectroscopically identified contaminants and non-globular cluster fillers. Columns 1 and 2 give the position in right 
ascension and declination (J2000), respectively. Columns 3, 4 and 5 are the Subaru/Suprime-Cam photometry. Measured recession 
velocities, when available, can be found in column 6. Column 7 lists the type of contaminant. The full table is available online. 
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Table |6] shows the position of those contaminants and their 
recession velocity. 

The spatial distribution of kinematically confirmed GCs 
is shown in Fig. [1] We look for rotation in the GC system 
around NGC 4494. We extend the kinemetry method de- 
scribed in Section 13.1.21 to discrete velocities and perform 
rolling or moving radial fits for the amplitude of the rotation 
(K-ot), the velocity dispersion (a) and the kinematic position 
angle (PAkin) simultaneously. Rolling fits are similar in prin- 
ciple to rolling or moving averages (,Sawilowskvll2007l) and 
were also used in jPOOl (see also iKissler-Patig fc Gebhardtl 
1 19981 . who use a very similar method). They are performed 
by first using the inner A^'bin GCs to fit the kinemetry, then 
removing the innermost point and adding the next further 
point to refit the kinemetry, and so on, until all GCs have 
been exhausted. We use bins of 25 and 20 for the blue and 
red GCs, respectively. Gaussian line-of-sight velocity distri- 
butions are assumed. In practice, for the j"^ radial bin we 
minimise the likelihood ratio (Aj): 



E 



{a' + (AKbs.O') 



+ ln(a| + (Al/obs,»)') 



(7) 



where Kiod,i,j is again given by Eq. [3]with ijkin.j = gphot = 
0.87. Symbols PAi, Vohs,i and AV^obs.z are the position angle. 



recession velocity and uncertainty on the recession velocity 
for the z"^ GC, respectively. We assume that the kinematic 
axis ratio of the GC system is equal to the photometric axis 
ratio of the galax;y light because the GC kinematic data do 
not constrain the kinematic axis (i.e., ?kin,j) ratio well. 

Uncertainties on the fits to the GC kinematics are ob- 
tained using a bootstrapping method similar to that used 
by Cote ct al. (2001). Basically, we obtain 1000 "mock" GC 
kinematic samples of 109 GCs (our actual sample size) by 
sampling with replacement from our measured distribution 
and fit each mock sample. The quoted uncertainties are the 
68% confidence interval inferred from the mock fits. 

We also notice that this method tends to enhance the 
estimated rotation value when the position angle varies as 
a free parameter. We used Monte Carlo methods to quan- 
tify this bias. We find that Vrot estimated using the above 
method on a sample of the same bin size and known slow 
rotation (Kot < 10 km s~^) can be enhanced by as much as 
20 km s^^ when PAkin is allowed to vary freely. This bias 
sharply reduces to < 5 km s~^ for known input Kot > 50 
km s~^. For this reason, fixing PAkin to a reasonable value 
is recommended whenever possible to avoid this bias, espe- 
cially when the output Kot is low. 

We find that the kinematic position angle for the blue 
GCs (170^14 deg) is consistent with the photometric po- 
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Figure 12. GC velocity distribution of the blue (hollow blue symbols), red (filled red symbols) GCs and UCDs (star symbols) recession 
velocities as a function of position angle (left panel) and galactocentric radius (right panel). Solid and dashed blue/red lines show the 
kinematic fit results and ±2(7, respectively, for the blue/red GC subpopulation. Blue/red dotted lines show 2x the velocity dispersion 
(corrected for rotation) intervals for the blue/red GC subpopulation. Black dashed line represents NGC 4494's systemic velocity. The 
photometric position angle is PA]^^^ = 173 deg. This figure is available in colour in the online version. 



sitio n angle in the inner parts of the galaxy measured by 
iNOgI (i.e., PAphot ~ 178 deg, see Fig. |6}. As will be dis- 
cussed in great detail below, the red GCs show lower rota- 
tion, making the determination of the kinematic PA highly 
uncertain and artificially enhancing Kot. Therefore, we fix 
the position angle to the photometric position angle (i.e., 
-P^kin.j = PApiiot) for both subpopulations in order to bet- 
ter constrain the other parameters. The resulting fits for our 
fiducial colour split are shown in Fig. [T2]for the whole sam- 
ple and rolling fits as a function of radius can be found in 
Fig. [6] for both GC subpopulations. We exclude GC88 from 
the radial rolling fits as its position at > 8re skews the mea- 
sured rolling radius significantly. This does not change the 
amplitude of the fitted parameters significantly. We find sig- 
nificant major-axis rotation at all radii for the metal-poor 
GCs only with 41 < Vrot < 95km s"^ (see Fig. ^. The 
metal-rich subpopulation may not show significant rotation 
as Vrot is consistent with km s~^ at most radii. To test the 
effect of our colour split assumption, we perform rolling fits 
as a function of {g — i)o colour with moving bins of size 25, 
which we show in Fig. 1131 These demonstrate that interme- 
diate colour GCs (i.e., 0.90 < {g~i)o < 0.98), in particular, 
do not show significant rotation and that this results in a 
higher uncertainty on the kinematic position angle. This pe- 
culiarity is yet another hint of a possible third intermediate 
colour GC subpopulation around NGC 4494. On the other 
hand, the reddest GCs {{g — i)o > 0.98) tentatively show 
some rotation. In every case we find that the velocity dis- 
persion (a) of all subpopulations are consistent with each 
other (see Fig. 



I 

E 100 
50 



> 



J 50 
MOO 
^ 50 


?300 



<200 - 



100 



'fixed' PA 




■ I I I I I I I I I I I I I I I I 




" i I I I I I I I I I I I I I I r 







vdria'ble' PA 




J I I I I I I I I I I I I I I I . 




I I I I I I I I I I I i i ii I I r 



,_L_ 



0.8 0.9 

(g-i)o 



0.8 0.9 1 

(g-i)o 



Figure 13. Rolling colour fits to the GC kinematics. Two cases 
are presented; PA]^[^ = P^phot (left panels) and variable PAj^jn 
(right panels). Shaded areas represent 68 per cent confidence in- 
tervals. 



3.2.4 GC metallicities 

We are able to measure spectroscopic metaUicities for 54 
individual GCs and the three UCDs aroun d NG C 4494 us- 
ing the CaT and the method described in iFlOl . Briefly, we 
fit the GC spectra using the pPXF routine described in 



ICappellari fc EmsellemI (|2004l ) to obtain a template of our 
GC spectra. We use the set of templates described in Sec- 
tion |323] (example fits are shown in Fig. [2)). Next, the fit- 
ted spectra are continuum normalised using the iraf rou- 
tine CONTINUUM with a splineS function of order ~ 4 and a 
stricter lower sigma clipping to av oid sp ectral features be- 
ing fitted. The CaT index defined in lFlol fCaTpin. see Table 
[3| is measured on the template fitted and continuum nor- 
malised spectra. We emphasise that this CaTpio is differ- 
ent from CaTpog used above for galaxy spectra, which has 
broader passbands to accommodate the broadening of the 
CaT features due to large velocity dispersion in galaxies. 
The CaTpio indices are transformed into metallicity using 
the following equation: 



[Fe/H] = -3.641 + 0.438 x CaTpio (8) 

based on the empirical conversion of lArmandroff fc ZinnI 
1 19881 ) derived from the integrated light spectra of Galac- 
tic GCs. For more information on this technique, see lFld . 

Fig. 1 141 shows our measured CaTpio index as a function 
of {g — i)o colour. We keep multiple CaTpio measurements 
of indivi dual obje cts as separate dat a poin ts. Predictions 
from the and iBruzualfc" CharlotI (|2003l ) simple stellar 
population models are overlaid for comparison. 

There are several immediately striking features in Fig. 
1141 First of all, the correlation between the CaTpio and 
colour is obvious and consistent with linear, albeit with 
large observational scatter. There is one outlier, namely the 
faint (io = 21.8) red GC, GC102. There is nothing obviously 
wrong with the photometry or spectrum (fitted or raw) of 
GC102, so we cannot explain the position of this (low signal- 
to-noise) outlier in Fig. 1141 

The bulk of th e GC data lie above the 
iBruzual Sz CharlotI (|2003|) single stellar population models 
and close to those oflV03! (for bluer colours). The appar- 
ently linear relationship found in NGC 4494 data u p to the 
reddest colours are in contrast to the findings of IFIOI for 
the giant elliptical galaxy N GC 1 407, where the generally 
redder data follow the 13 Gyr Ivo^ model track more closely 
at all colours. Many blue/red GCs in NGC 1407 have 
higher/lower measured CaTpio indices than the NGC 4494 
GCs. We discuss this discrepancy and its implications for 
using the CaT as a metallicity indicator for extragalactic 
GCs in Appendix IbI 

Finally, there appears to be a concentration of GCs with 
CaTpio ~ 6.5A (or [Fe/H]~ -0.9). While our sample of 
spectroscopic GC metallicities is rather small for inferring 
the global properties of the distribution of GC metallicities 
in NGC 4494, it appears to be single peaked despite exhibit- 
ing clear {g — i)o colour bimodality (KMM yields a p- value of 
0.002 for the confirmed GC sample). We not e however that 
KMM is less reliable for small samples (see lAshman et al.l 
1 19941 ). Nevertheless, it is puzzling that the peak in the CaT 
distribution corresponds to the trough of the colour distribu- 
tion (see Fig. 1141). Th is res ult is reminisce n t of th at obtained 
for NGC 1407 (|F10|) and ICaldwell et al.1 (j201ll ). where the 
clear colour bimodality also translated into a skewed single- 
peaked spectroscopic metallicity distribution. 
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Figure 14. CaTpio and inferred metallicity (right y-axis) as a 
function of colour for GCs (filled circles) and UCDs (hollow tri- 
angles) around NGC 4494. All repeat CaTpio measurements of 
individual objects are shown. Histograms for each axis are also 
shown. One outlier (highlighted with a large hollow circle) is dis- 
cussed in the te xt. Predicti ons from the single s t ellar population 
models of IVOSl (black) and iBruzual fc Charloti ||2003| . grey) are 
shown for 5 (dotted line), 9 (dashed line) and 13 (solid line) Gyr. 



3.2.5 Three UCDs around NGC 4494 

As briefly mentioned in Section 13.2.11 we report 

the disc overy of 3 spectroscopically confirmed UCDs 

( Drinkwa ter et al.l bood ) associated with NGC 4494. We 
adopt a magnitude definition for UCDs of Mj < —10.9, 
roughly equivalent to that adopted by lEvstigneeva et all 
(200i) in the V^-band. The 3 UCDs have absolute i-band 
magnitudes of -11.88, -11.92 and -11.71 mag, well within 
the range for UCD luminosities. Most UCDs have typically 
been found in dense cluster environments although, for ex- 
ample, one UCD has been confirmed around the Sombrero 
galaxy (|Spitler et al.ll2006l : lHau" et al. 1(2009), a spiral galaxy 
in a low-density environment. Another relevant example is 
that of NGC 5128, also an L* early-type gal axy, wherein 
sever al possible UCDs may have been found (|Tavlor et al.l 
[2O1O). 

There is ongoing discussion in the literature about the 
origin and definition of UCDs. Popular formation scenar- 
ios propose that UCDs are either (1) the bright end of the 
compact star clus ter (i.e., GC) luminosity function (e.g., 
iMie ske et a"l]|2004l ). or (2) the slightly more extended rem- 
nants of tidally stripped dwarf galaxies (Bekki o t al. 20 03|) . 
Of course, both scenarios may oc cur (see iDa Rocha et al.l 
I2011I : n orris fc KannappanI 1201 ll ). We thus examine the 
properties of the UCDs around NGC 4494 in order to deter- 
mine their most likely origin. 

UCDl and UCD2 are within the HS T imaging. We 
measure their sizes using ishape (iLarsenl [l999l ) and ob- 
tain 1.7 ± 0.4 and 2.3 ± 0.6pc for UCDl and UCD2, re- 
spectively. We are unable to obtain a size estimate for 
UCD3 as it is unresolved on the Subaru image. There- 
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fore, at least two of the three UCDs around NGC 4494 
are compact, as predicted by the star cluster origin sce- 
nario. All three UCDs have high CaT-inferred metallic- 
ities (—0.3 > [Fe/H] > —0.4) as is also th e case for 
some other UCDs (e.g., lEvstigneeva et aLlbOOTl). UCDl is 
detec ted with Chandra in the X-ray ( Humphrey fc Buotd 
|2008| . object id NGC_4494_CXOU_12:3125. 5-^254619) with 
Lx = (2.1 ± 1.1) X 10^* ergs s"^ This X-ray luminosity is 
consistent with the presence of low mass X-ray binary stars 

in UCDl. 

FoUowing lDa Rocha et al] (|201ll ). we compute the num- 
ber of expected GCs brighter than io = 19.3 (i.e., the mag- 
nitude of our faintest UCD candidate) around NGC 4494 
assuming that UCDs around NGC 4494 are simply the 
bright end of its GC system. We use the total number of 
GCs and the GC luminosity function (|Kundu fc Whitmord 
[iool) to estimate that 2 ± 2 UCDs are consistent with be- 
ing the bright extension of the GC system. Thus, under 
these assumptions, our 3 confirmed UCDs brighter than 
io = 19.3 around NGC 4494 are consistent with repre- 
senting the bright end of the GC luminosity function. We 
conclude that there is no need to invoke another forma- 
tion channel such as t i dal stripping of dwarf galax ies (e.g., 
[dT Rocha et al. I I2OIII : iNorris fc KannappanI I2OI0I ) for the 
UCDs around NGC 4494. 



4 DISCUSSION 

In this section, we compare the observed properties of 
NGC 4494 with the predictions from theoretical models in 
order to get an understanding of its formation. 



4.1 Inferences from the stellar light 



In Section 13.1.21 we report a 'flattening' of the stellar kine- 
matics of NGC 4494 with radius such that the kinematics 
become more disky at large radii. In other words, at large 
radii only the stars close to the semi-major axis show ro- 
tation such that the kinematic axis ratio ^kin is low. This 
is supported by the observed anti-correlation between /13 
and (V^bs — VsyB)/^, which indicates the presence of a disk- 
like structure. However, it contrasts with the stellar surface 
brightness of the galaxy, which has a very constant and rel- 
atively round profile at all radii. This kinematic fiattening 
at large radii may be related to the transition suggested by 
iHoffman etall (I2OI0I ) between 1-3 r^, where the kinematic 
signature of the progenitors' disk stars survived. However, 
it is a puzzle as to how such a 'kinematic' disk at large radii 
could be invisible in the imaging data. Indeed, the axis ratio 
of kinematic sub-structures are usu ally found to agree wit h 
that of the stellar distribution fe.g.. iKrainovic et al.|[200^ . 

The stellar populatio ns don't show a ny hi nt of a fiat- 
tened distribution either. iDenicolo et al.l l|2005l ) reported a 
central age of 6.7 Gyr with a central metallicity of [Fe/H] ~ 
-1-0.03 dex for NGC 4494 and Hoi find no evidence for ra- 
dial metallicity variations with metallicities roughly con- 
stant around [Fe/H] > —0.5 dex between ~ 0.2 and lAr^. 
In this work, we find no azimuthal variations in the mea- 
sured CaT index values (see Fig. [S| or colour. The CaT 
gradient suggests that the luminosity weighted metallicity 
of the NGC 4494 stars is either higher than ~ —0.5 dex 



from 0.2re all the way out to at least ~ 1.6 and/or that 
the metallicity gradient is undetectably shallow for the CaT 
method. The colour gradient in the same radial range sug- 
gests a moderate metallicity gradient of ~ —0.17 ±0.02 dex 
per dex. 

Recently, the SAURON team reported the results of 
thei r 2D stellar population analysis on 48 early-type galax- 
ies (|Kuntschner et al]|2010l ). They fin d that flattened struc- 
tures in the images o f fast-rotators (jEmsellem et al.ll2007l : 
IKrainovic et al.ll2008h with disky kinematics have distinct 
stellar populations, while galaxies classifled as slow rota- 
tors, and sometimes harboring inner kmematically decoupled 
cores, show no clearly distinct stellar population variation. 
NGC 4494 does not appear to fit either of those categories. 
At small radii ~ O.lre , the kinematically decoupled core 
suggests it is a slow-rotator. At all other radii, its increas- 
ingly fiattened kinematics suggest it is a fast rotator, but 
as stated above we find no evidence for a flattening of the 
stellar distribution or of distinct stellar populations from ei- 
ther the CaTpo9 or the modelling of the surface brightness 
profile and colours. In any case, transitions between slow 
and fast rotators at large radii may be common as £09. also 
report a transition from a fast to a slow rotator beyond the 
SAURON field-of-view in NGC 821. Again, these may be the 
first observational evidences for the transitional kinematics 
expect ed to occur between 1 and 3re in major merger rem- 
nants (jHoffman et al.ll2010l ). 

In Fig. [TSl we plot a standard (Kot/c r) vs e = 1 — gphot 
diagnostic diagram (jCappellari et al ]|2007l ). This plot can be 
used to diagnose both the intrinsic structure of the galaxy 
(in particular NGC 4494) and the nature of the merger. 
We have highlighte d NGC 4494 on this diagram showing 
the ATLAS3D data l|Emsellem et al.ll2011^ . The green curve 
in Fig. [TJ] represents an oblate isotropic rotator seen edge- 
on, and the magenta curve shows what is typical for an 
edge-on fast rotator after mo delling the dynamics o f the 
SAURON sample as inferred bv lCappellari et al.1 l|2007l ). Be- 
cause galaxies are generally observed below the green curve, 
this suggests that they are not isotropic. These results al- 
low us to derive a best-guess solution for the inclination of 
any fast rotator (including NGC 4494) by assuming that 
its internal anisotropy follows the mean trend of the other 
galaxies. We can also thereby estimate the ellipticity and 
(Kot/c) values that would be obtained if the galaxy was 
viewed edge-on. 

The black curve shows the track of possible edge-on val- 
ues for NGC 4494 for a series of different inclinations, where 
we note that the dependence of the dispersion on inclination 
enters through the assumed anisotropy. The intersection of 
the black curve with the magenta curve then represents the 
self-consistent solution for NGC 4494 under the SAURON- 
based anisotropy assumption. From this we conclude that 
NGC 4494 is most Itkely (but not definitively) a flattened 
galaxy (Qphot ~ 0.6) seen at an inclination of ~ 45 degrees. 
It may even be an SO rather than a bona fide elliptical. 
The uncertainties here are driven not by the measured pa- 
rameters (which are determined very precisely) but by the 
intrinsic scatter in the anisotropy-ellipticity relation. To es- 
timate this, we use the spread of observed galaxies to the 
right of the magenta curve in Fig. 1151 This suggests an in- 
trinsic ellipticity uncertainty of ~ ±0.2, and an inclination 
between 39 and 90 degrees. Similarly, using 2D dynamical 
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Figure 15. Azimuthally averaged rotation (V^o t/c-) versus ellip- 
ticity (e = 1 — gphot) diagnostic diagram, after ICappellari et alj 
||2007| ) .Data points represent the central regions of nearby early- 
type galaxies classified as "fast-rotators' from the SAURON and 
ATLAS3D surveys. The green dashed curve shows a theoreti- 
cal prediction for edge-on oblate isotropic rotators, and the ma- 
genta dotted curve shows the inferred edge-on average relation 
for nearby fast-rotators. The observed position of NGC 4494 in 
the diagram is labelled, with a solid black curve showing the path 
of possible intrinsic values for a sequence of different assumed in- 
clination angles. This figure is available in colour in the online 



modelling of NGC 4494, iRodionov fc Athanassoulal (|201ll ) 
found that a ~ 45 degrees inclination may be preferred. If 
this interpretation is correct, then to recover edge-on Vrot 
estimates, all of our velocity estimates should be increased 
by ~ 40 per cent. 

4.2 Inferences from the GC system 

We find that the GC colour subpopulations are reasonably 
well defined in NGC 4494 and choose a nominal colour cut 
at (g — i)o — 0.99 to delineate metal-poor from -rich GCs. 
The 54 spectroscopically measured GC metallicities vary be- 
tween -2.0 dex < [Fe/H] < 0.0 dex and their distribution 
appears single-peaked around [f e/Z/J^ean ~ —1.0 dex de- 
spite their clear bimodal colour distribution. A s imila r be- 
havio ur was found for GC s around NGC 1407 (|f10| ) and 
M31 (|Caldwell et al.|[201ll ). These may have some interest- 
ing implications for the ubiquity of the GC metallicity bi- 
modality as inferred from GC colour distributions (also see 
lYoon et al. I I2OO6I : iPeng et al.ll2006l : iBlakeslee et aLlbOld l. 

It has been suggested (and shown) that the distribu- 
tion of red GCs usually follows that of the galaxy stars, 
while blue G Cs follow the X-ray halo profile of their host 
galax y (e.g., iMinnitil \l99^ : iForbes et all |2004| : iBolev et all 
l2009h . This suggestion is based on both model predictions 
and the similarities of the respective spatial distribution and 
abundances. Surprisingly, we find that the surface brightness 



profile of the stars does not compare well with the spatial 
distribution of the red (or the blue) GCs in NGC 4494 (see 
Fig. Ilip . Indeed, the Sersi'c index of the stars is inconsistent 
with that of both GC subpopulations, and this difference is 
even more pronounced for the red GCs. A similar conclusion 
is reached from comparing the kinematics of the stars, which 
agrees better with the blue GC kinematics than with that 
of the red GCs at the same radius (see Fig. [G]). On the other 
hand, the colour of the galaxy stars is consistent with that 
of the red GCs (see Fig. [8]). Therefore, the association of 
the red GCs with the galaxy stars is less clear in NGC 4494 
than previously observed in other galaxies. 



In Section 13.2.31 we measure the rotational velocity 
and velocity dispersion of the blue and red GCs. In or- 
der to put our results into a broa der contex t . Fig. 1161 re- 
produces parts of figure 12 from iLee et all (|2010l ) usmg 
the updated kinematic data presented in Tableland com- 
pares GC kinematics with the global properties of the host 
galaxies for large ETGs. One caveat of the following com- 
parison has to do with the heterogeneity of the methods 
used in the various dynamical studies represented in Fig. 
1161 Moreover, NGC 4494 is the only fast-rotator with large 
GC kinematic sample. Thus, its GC kinematics may be in- 
trinsically different. Nevertheless, we find that NGC 4494 
compares well with other large ETGs. We confirm that it 
agrees with the trends found between the velocity disper- 
sion of the whole GC system and X-ray, central galactic 
velocity dispersion and j3 -band absolute magnitude as re- 
ported in lLee et al.l l|2010l ). However, the kinematics of the 
blue GCs of NGC 4494 appear to deviate from that of 
other large elliptical galaxies in the Frot.biuc/crbiuc vs astar 
space, emphasising the large rotation of the blue GCs. Per- 
haps surprisingly given the arguably "unusually low" X- 
ray luminosity of NGC 4494, the GC kinematics agree well 
with the expected trend with galaxy X-ray luminosity (e.g., 
iRomanowskv et al.|[200l : IO'Sullivan fc Ponmanll2004i v This 
also suggests that the processes involved in the formation 
of NGC 4494 that led to its peculiar X-ray luminosity and 
kinematics have preserved the majority of the correlations 
between its global properties and its GC kinematics as with 
other large ETGs. 

Assuming that NGC 449 4 has indeed undergon e a re - 
cent interaction as inferred bv lO'SuUivan fc Ponmat] (|2004l ) 
from its X-ray luminosity, then we may reasonably conclude 
that the bulk of its GC system formed before that inter- 
action and is thus a conglomerate of the progenitors' GC 
systems. Therefore, the current GC kinematics may hold 
clues to the understanding of the kinematics of the progeni- 
tor galaxies involved in the latest merger event. In this spirit, 
we use numerical simulations of disk-disk major mergers in 
order to see what can be learned about the progenitors from 
the GC kinematics. Details of the simulations and results 
are presented in Appendix [C] By studying numerous simu- 
lations of mergers with a variety of initial conditions, we find 
that the model that best reproduces the observed kinematics 
of the blue and red GCs (ignoring the intermediate colour 
GCs that show little rotation) involves a major disk-disk 
merger with a large amount of orbital angular momentum 
and a retrograde-retrograde orbital configuration. Models 
with other orbital configurations do not yield GC kinemat- 
ics that are consistent with what is observed for NGC 4494. 
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Figure 16. Comparison of salient GC kinematic properties with 
host galaxy global parameters: X-ray luminosity (log(Lx))i stel- 
lar velocity dispersion (a star) and absolute V-band magnitude 
{Mv)- Symbols on the y-axis and GC kinematics are from Table 
[5l We include literature GC kinematic studies (filled circles) and 
NGC 4494 data (hollow stars and upper limit sy mbols). Sol i d lines 
are bisector fits to correlated parameters from I Lee et akl 1I2OIC1I . 
figure 12). Datum for NGC 4494 disagrees with the Vpot blue/^blue 
vs Ostar trend observed in other large elliptical galaxies (middle 
lower panel). 



This is a suggestion that NGC 4494 may be a major merger 
remnant. 



5 SUMMARY AND CONCLUSIONS 

NGC 4494 has often been dubbed an 'ordinary elliptical' 
galajcy, making it an ideal target for understanding the for- 
mation and evolution of a typical galaxy. In this work, we 
combine imaging and spectroscopy of stellar light and GCs 
in NGC 4494. From the imaging we obtain the stellar surface 
brightness profiles and GC density profiles. We find that the 
colour distribution of the GC system around NGC 4494 is 
statistically better fitted by two Gaussians than a single one 
suggesting metallicity bimodality. The spectroscopy yields 
spatially resolved kinematics and abundances of the stellar 
component and GCs. A total of 109 individual GCs and 3 
UCD/DGTOs are confirmed spectroscopically. The proper- 
ties of the UCDs are consistent with them being bright GCs. 
Metallicities are measured for 54 GCs and the 3 UCDs us- 
ing the CaT absorption lines. The CaT inferred metallicity 
distribution for the GC system is single-peaked, in contrast 
to that inferred from the colours, which is clearly bimodal. 
Measurements for both the stellar light and GC spectroscopy 
can be found in Appendix [X] We find that while the inter- 
mediate colour (green) GCs do not rotate, the blue (and 
possibly reddest) GCs do rotate (Fig. I13p as do the galaxy 
stars (Fig. |6]). The velocity dispersion of all GC subpopula- 
tions and the galaxy stars are consistent. 

A comparison of the GC system's kinematics and the 
global properties of NGC 4494 with those of other large 
ETGs suggest that NGC 4494 is indeed typical with the 
exception of the unusually high blue GC rotation. We find 
suggestive evidence in the distinct kinematics for a possible 
third intermediate colour (green) GC subpopulation. We in- 



fer that most of the observational evidence suggests that 
NGC 4494 is consistent with formation via a recent gas-rich 
major merger, but other formation scenarios cannot be ruled 
out. Some remaining open issues are: 

(i) the relatively old (6.7 Gyr) measured central age sug- 
gests that little dissipation and recent central star formation 
has occurred (i.e., possibly a dry merger). However, its cuspy 
surface brightness profile and the presence of a kinematically 
decoupled core suggest a high gas fraction of > 15 per cent 
(Hopkiiis ct al. 200 9; Hoffm an et al., 2Q1Q ). This is also sup - 
ported by the presence of an inner dust ring (|Xu et al.ll2010l ) . 
Therefore, the evidence is not conclusive as to the amount of 
gas dissipation involved in the formation of NGC 4494. The 
possibly 'unusually' low X-ray luminosity suggests that the 
mechanism involved in depleting the gas must have allowed 
enough dissipation to form a kinematically decoupled core 
without significant recent central star formation. 

(ii) The inner kinematic profile of NGC 4494 follows the 
light profile well up to ~ 1.8re, beyond that the kinematics 
become increasingly flattened with radius (i.e., the axis ratio 
gkin is low). We find no evidence of a flattened component in 
the stellar distribution from our imaging. It is however not 
surprising to find transitions in the kinematics of merger 
remnants at l - 3re a s they are predicted in the models of 
iHoffman etHI l|2010l ). We suggest that NGC 4494 may be a 
flattened galaxy (qphot = 0.6±0.2), possibly even an SO, seen 
at an inclination of 45lg^ degrees based on its position in the 
classic (Kot/o") vs e = 1 — Qphot diagram (Fig llSp . This may 
explain its high flattened rotation with round photometric 
isophotes. 

Off-axis deep optical spectroscopy out to large galac- 
tocentric radii would be a good way to independently test 
the ff aliened kinematics. It would also allow for a stellar 
population analysis (i.e., ages, metallicities and a-element 
abundanc es) using the standard a nd well tested Lick sys- 
tem (e.g.. IProclor fc SansomI [20021 ). The added knowledge 
of spatially resolved stellar ages and a-element abundances 
may highli ght distinct stellar pop ulations as observed by 
SAURON (jKuntschner et al.ll2010l ) that were not seen from 
our near-infrared spectra. 

This work shows that such complete studies of individ- 
ual galaxies incorporating as many lines of evidence as avail- 
able can help disentangle the processes involved in the for- 
mation and evolution of selected galaxies. As such datasets 
become available, observational constraints based on statis- 
tical samples will help improve our understanding of the 
details of galaxy formation and evolution. 
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APPENDIX A: DATA TABLES 

This section contains our measurements for the photome- 
try, kinematics and stellar populations of stars and GCs in 
NGC 4494. 

APPENDIX B: THE CAT AS A METALLICITY 
INDICATOR 

In this section, we focus on the reliability of the CaT as a 
metallicity indicator for unresolved populations (i.e., GCs 
and galaxies). We use the data presented in this work to 
shed some lig ht on two open issues. 

First, the lV03l SSP models predict that the CaT satu- 
rates at high metallicity (i.e., [Fe/H\ — —0.5). This predic- 
tion is consist ent w ith the distribution of the GC data for 
NGC 1407 in iFld but is less clear from Fig. [14] where the 
relationship between colour and CaTpio appears linear at all 
probed metallicities. This is con sistent with what has bee n 
observed in the Galactic GCs (|Armandroff fc ZinnI Il98^ \ 
On the other hand, in Fig.[8]the measured CoTpog for galaxy 
spectra scatters about a constant val ue, wh ich coincides with 
the saturation value predicted by the lV03l models. This may 
be interpreted as proof of the saturation prediction. How- 
ever, the galaxy colour gradient (see Fig. \^ is shallow at 
similar radii suggesting that the expected change in CaTpog 
may be small. In other words, the absence of a Ca T radial 
gradient i n the galaxy found in this work and in F09 may not 
prove the lVOSi saturation prediction. The question remains 
open. 

Second, iFld found that the brightest blue and red GCs 
in NGC 1407 have the same CaT index value despite their 
wide separation in colour. This cast serious doubt on the 
reliability of the CaT as a metallicity indicator. Moreover, 
several fitted GC spectra across the range of GC luminosities 
probed around NGC 1407 showed Paschen lines. We find no 
such Paschen lines in the GC spectra for NGC 4494, which 
cover a comparable range in absolute luminosities. The pres- 
ence of Paschen features in the GCs around NGC 1407 can- 
not be confirmed directly on the raw spectra. 

Several possible explanations for these behaviours are 
put forward in FIO, including the presence of hot blue stars 
such as blue horizontal branch, blue straggler or young stars 
mainly in the blue GCs whose Paschen lines might be affect- 
ing the CaT features. Another possibility is that the CaT 
may saturate at lo wer m etallicity (i.e. [Fe/H] ~ —0.8 dex) 
than predicted by Ivosl or perhaps colour does not trace 
metallicity linearly (e .g., lYoon et al.ll2006l : iPeng et al.ll2006l : 
iBlakeslee et~alll2010l ) . 

Based on the results for NGC 4494 shown in Fig. [HI 
there is no clear evidence that the CaT saturates earlier than 
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Table Al. Individual values for NGC 4494 galaxy light. Columns 1 and 2 give the position of the individual slits in right ascension and 
declination (J2000), respectively. The observed veloc ity moments Vobsi and hn, appear in columns 3, 4, 5 and 6. Measured values of 

the CaT index and [Fe/H] from the method of IF09| . when available, are shown in columns 7 and 8. We quot e the maximum metallicity 
(i.e., [Fe/H] = +0.2) when the measured CaTpog is higher than the metallicity range available in the |V03| models. The errors on the 
quoted [Fe/H\ values are based on the index errors only and do not include possible systematics due to the adopted calibration. The 
full version of this tabic is available through MNRAS. 
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Table A2. Individual spectroscopically confirmed GCs and UCDs around NGC 4494. Columns 2 and 3 give the position in right ascension 
and declination (J2000), respectively. Columns 4 to 10 present the photometry. The measured recession velocity is given in column 11, 
while values of the CaTpio index and inferred [Fe/H], when available, are in columns 12 and 13, respectively. The errors on the quoted 
[Fe/H\ are based on the index errors only and do not include systematics e.g. adopted calibration. The full version of this table is 
available through MNRAS. 
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predicted bylVOd or that the CaT behaves non-linear ly with 
colour. In fact, the relationship between {g — i)o colours and 
CaTpio is consistent with being linear, albeit with large ob- 
servational scatter. This suggests that the strange distribu- 
tion of CaTpio values in the GCs around NGC 1407 may be 
best explained by the presence of hot blue stars in a signifi- 
cant number of its GCs afterall. 



APPENDIX C: NUMERICAL SIMULATIONS 
OF THE GC SYSTEM 

In this section, we discuss whether the observed kinematics 
of blue and red GCs around NGC 4494 can be reproduced 
reasonably well via major merging between two disk galaxies 
with pre-existing blue and red GCs. 



CI Model description 

In the following, we ignore the presence of the non-rotating 
intermediate colour GCs. We run dissipationless simulations 
(i.e., no gas dynamics and no new formation of stars and 
clusters) and compare the physical properties of the GC sys- 
tems in the merger remnants to search for the model that 
best reproduces the following observed properties of blue 
and red GCs: 

(i) the maximum rotational velocity ('Kot,max) is as large 
as 100 km s"'^ for blue GCs, 

(ii) the central (maximum) velocity dispersion (a) is as 
large as 150 km s^^ for both subpopulations, 

(iii) Vrot,max/o" is larger than 0.3 for both metal-poor and 
-rich GCs, and 

(iv) Kot.max is larger in the blue than in red GCs. 

We consider that the latest (iv) is one of the key physical 
characteristics of the GCs around NGC 4494. 

Since the numerical methods and techniques we em- 
ploy for modelling the dynamical evolution of dissipationless 
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Figure CI. Line-of-sight rotational velocities (Vrot, solid lines) 
and velocity dispersion {a, dotted lines) for blue and red GCs 
along the x-axis in the merger remnant for the fiducial model 
(Model 1). Here the y-components of GC velocities are used for 
deriving the line-of-sight Vrot and a at each bin along the x-axis. 
This figure is available in colour in the online version. 

mergers between two disks with GCs h ave been detailed else- 
where (.Bekki et al. 2005 ; Bekkill2010l ). we give only a brief 
review focusing on the main particularities of this project 
here. The progenitor disk galaxies taking part in the merger 
are given a dark matter halo, bulge, thin exponential disk, as 
well as GCs initially located in the bulge (hereafter BGCs) 
and in the halo (HGCs). The total mass and size of the 
exponential disk (bulge) are A/d (Afb) and _Rd {Rh), respec- 
tively. We show results from models where the progenitor 
disk galaxies are similar to the Galaxy (i.e., Afb = 0.167Afd 
and i?b = 0.27?d). In this work, the total (virial) mass of the 
dark matter halo (Afdm) in a disk galaxy is set to be 9Afd. 

The initial distribution of HGCs in the progenitor disk 
galaxies follows a power-law profile with index a = —3.5. 
The extent of the HGC system (-Rhgc) is assumed to be 
3-Rd. In other words, the power-law distribution is truncated 
at _R = 3i?d. The half- number radius of the HGC system 
is set to be 0.29-Rd. The initial distribution of BGCs also 
follows the same power-law but has a truncation radius of 
fgc^^HGC, where rgc controls the compactness of the spa- 
tial distribution of BGCs. We find that rgc is important for 
determining the kinematical differences between HGCs and 
BGCs in the merger remnants. 

In the present models, the inner metal-rich BGCs and 
outer metal-poor HGCs ultimately become blue and red 
GCs in NGC 4494, respectively. This is a simplification since 
additional metal-rich GCs may form from gas dynamics dur- 
ing galaxy merging (e.g., iBekki et al.ll200^ ). Even so, the 
models may enable us to understand how large amounts 
of rotation are present in metal-poor and possibly -rich 
GCs around NGC 4494. For convenience, K-ot,max (cr) of 
HGCs (blue GCs) and BGCs (red GCs) are referred to as 
Kot.max.biue (cTbiue) and K-ot,max,red (cTred), respectively. Here 
Rd and Md are set to be 13.4 kpc and 4.3 x lO'^'^M©, respec- 
tively, to reproduce the observed kinematics of the stars in 
NGC 4494. 

The mass ratio of the two progenitor disks in the sim- 
ulated mergers is 1 (i.e., equal- mass major merger) and the 
orbit of the two disks is initially set to be in the a;j/-plane 
for all models. The pericenter distance (rp) and the eccen- 
tricity (ep) of the merger are free parameters that infiuence 
the orbital energy and angular momentum. The spin of each 



Figure C2. Same as Fig. ICll for the radial dependence directly 
comparable with the observational results in this work. Vertical 
dotted line indicates the effective radius of NGC 4494. This figure 
is available in colour in the online version. 

progenitor galaxy is specified by two angles 6i and 4>i. 9i is 
the angle between the z-axis and the angular momentum 
vector of the i"" progenitor disk. is the azimuthal angle 
measured from the x-axis to the projection of the angular 
momentum vector of the i*** progenitor disk onto the xy- 
plane. 

We mainly show the results for our "fiducial model", 
which has rp = 2Ri, ep = 0.72, Oi = 210°, 62 = 225°, 
(pi = 45° , and 4)2 ~ 120° , as it best reproduces the observed 
four key properties of the NGC 4494 GCs listed above. 
The fiducial model has an orbital configuration similar to 
"retrograde-retrograde" merging (i.e., the orbital spin axis 
of the merger is anti-parallel to the intrinsic spin axes of 
the two disks). We also show the results of the models with 
"prograde-prograde" ( "prograde-retrograde" ) orbital config- 
urations in which 61 — 30° and O2 — 45° (^1 = 30° and 
02 = 225°) and all other parameters unchanged. Model pa- 
rameters and some salient results are given in the Table ICll 

When estimating the kinematics of the GC system and 
binned major-axis profiles, the merger remnant is assumed 
to be viewed near to edge-on. In order to have enough ob- 
jects in each bin, GCs at corresponding minor axis distances 
are included. 



C2 Results 

Fig. lCll shows the rotational velocity (Kot) and velocity dis- 
persion (a) along the x-axis of the merger remnant for the 
fiducial model. The y-component of GC velocities for each 
of 10 bins are used for deriving line-of-sight profiles. Fig. ICll 
reveals that both metal-poor and -rich GCs exhibit global 
rotation (—12 kpc < a; < 12 kpc) albeit with stochastic scat- 
ter due to the small number of GCs in each bin. In fact, the 
apparently rapid change in Vrot around x = 7 kpc is due 
largely to small number statistics. 

In order to compare the simulated profiles with the ob- 
served one, we investigate radial (R) profiles of the rota- 
tional velocity and velocity dispersion (Fig. IC2|I . In the fidu- 
cial model, the larger amount of orbital angular momentum 
(ep = 0.72 and rp = 27id) allows the outer components (dark 
matter halo and HGCs) to acquire a large amount of rotation 
due to the efficient conversion of orbital angular momentum 
into spin angular momentum. For R > Rc (=3.75kpc), the 
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Table CI. Model (column 1) parameter values and brief summary of results. Column 2 lists the pericenter distance of the merger in 
units of the disk size R^. The orbital eccentricity of the merger is shown in column 3. Description of the orbits in column 4 are coded as 
follows: "PP" , "PR" , and "RR" for prograde-prograde, prograde-retrograde, and retrograde-retrograde merging, respectively. Column 5 
and 6 shows the maximum rotational velocity of metal-poor and -rich GCs around the remnant, respectively. Columns 7 and 8 list the 
ratio of the maximum rotational velocity to the central velocity dispersion for metal-poor and -rich GCs, respectively. 



simulated Kot.biue can be larger than 50 km s ^ (see Fig. 
IC2|I in agreement with the observations presented herein. 

As shown in Table ICll the models with parabolic en- 
counters (i.e., models 2 and 3) that initially have a sig- 
nificantly smaller amount of orbital angular momentum do 
not show large Kot, max, blue compared to the fiducial model. 
This implies that the observed large rotational velocity 
(Kot/o" ~ 0.7) may help constrain the orbits of the pro- 
genitors of NGC 4494. Moreover, the small Cp in the fiducial 
model suggests that the progenitor galaxies were either a 
binary pair of galaxies or two large galaxies dominating a 
small bound group. 

Fig. IC2I also shows that even though the red GCs 
(BGCs) do have rotation, their rotational amplitude is 
smaller than that of the blue GCs. In the present models, 
the final kinematics of red GCs may become comparable to 
those of the stellar components of the merger remnants as 
a result of the adopted rgc = 0.1. The fiducial model has 
a retrograde-retrograde orbital configuration yielding a fi- 
nal stellar remnant with apparently low rotation (< 50 km 
s~^) depending on the viewing angle. Fig. IC2I shows that 
Kot,max,rod is siguificantly smaller than Kot,max,biue as ob- 
served in NGC 4494. Moreover, as shown in Table [CT] mod- 
els with prograde-prograde and prograde-retrograde orbital 
configurations do not yield large kinematic differences be- 
tween blue and red GCs. This suggests that the orbital con- 
figuration of disk-disk major mergers may be the key factor 
in causing the observed kinematic differences between blue 
and red GCs. 

We thus conclude that the observed kinematics of the 
GC system around NGC 4494 are broadly consistent with a 
formation scenario wherein the remnant galaxy formed via 
major disk-disk merging, involved large amounts of orbital 
angular momentum and an initial retrograde-retrograde or- 
bital configuration. An interesting prediction of the fiducial 
model is that the dark matter halo of NGC 4494 also has a 
significant amount of global rotation. The predicted maxi- 
mum line-of-sight rotational velocity of the dark matter halo 
is ^ 40 km s~^ and Kot/o" ~ 0.25. However, the simulated 
rotation profile of the red GCs (~ 50 km s"'^) is larger than 
the observed one (~ 20 km s~^) for ~ 1 — 2Rc (corre- 
sponding to 3.8-7.6 kpc). We do not find models that show 
comparably low remnant red GC rotation whenever the rem- 
nant rotation of the blue GCs is large (Kot, max, blue ~ 100 
km s~^). This could mean that physical processes other than 



the coUisionless major merger scenario explored in this study 
may be necessary to explain all observed stellar and GC 
kinematics NGC 4494 in a fully self-consistent manner. 



